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I.    INTRODUCTION 

Modern  weapons  systems  are  becoming  more  and  more  complex  and  are 
constantly  pushing  the  limits  of  technology.  Performance,  in  terms  of  agility  and 
faster  response,  are  demanded  at  the  same  time  that  reliability  and  operation  in 
a  wider  variety  of  environments  and  applications  are  required.  All  of  this 
highlights  the  importance  of  the  design  of  the  underlying  control  systems  of  the 
weapon. 

Many  times  a  choice  must  be  made  as  to  the  minimum  level  of  reliability  or 
robustness  acceptable  at  the  expense  of  a  certain  amount  of  performance.   For 
a  single-input  single-output  (SISO)  linear  quadratic  optimal  control  system  in 
which  there  is  no  cross-coupling  of  the  states  of  the  system  and  the  control 
inputs,  there  are  well  defined  and  guaranteed  levels  of  robustness.   However, 
for  a  multivariable  or  multi-input  multi-output  (MIMO)  system,  while  the 
guarantees  still  apply,  the  measurement  of  the  robustness  of  a  system  is  not 
necessarily  well  defined.  In  addition,  for  a  system  with  cross-coupling  of  the 
states  and  inputs,  which  can  occur  when  the  performance  measure  is  in  the 
generalized  quadratic  form,  there  are  no  longer  any  guarantees  on  the 
robustness.  An  optimized  system  of  this  type  may  be  vulnerable  to  even  minor 
uncertainties  or  deviations  from  the  model  used  to  design  the  control  system. 

Design  techniques  do  exist  for  working  with  these  problems,  but  they  are 
typically  of  a  graphical  nature  and  do  not  easily  provide  a  way  to  measure  and 
vary  the  trade-off  between  performance  and  robustness.   The  goal  of  this  thesis 
then  is  to  develop  a  design  method  for  using  one  variable  parameter  to  quickly 
and  easily  get  a  measure  of  this  trade-off.  The  focus  will  be  confined  to 


generalized  linear  quadratic,  finite-dimensional,  multi-input  multi-output  systems. 
It  is  assumed  that  full  state  feedback  is  available  so  observers  will  not  be 
necessary.  In  addition,  measurement  and  plant  noise  will  be  ignored  for  this 
development,  so  techniques  such  as  Kalman  filtering  will  not  be  necessary. 
After  developing  several  methods  for  designing  suboptimal  robust  controllers, 
these  will  be  simulated  using  a  statically  unstable  fighter  aircraft  which  provides 
a  convenient  MIMO  model.   Also,  the  required  performance  will  be  chosen 
such  that  the  cross-coupling  of  the  states  and  the  control  inputs  will  occur. 


II.  BACKGROUND  LINEAR  QUADRATIC   THEORY 

Standard  linear  quadratic  (LQ)  control  theory  is  well  understood  and  many 
tools  exist  for  designing  these  systems.  As  long  as  the  states  and  the  inputs  are 
not  coupled  in  the  performance  function,  these  systems  have  very  good  and 
well  defined  properties  of  stability  and  robustness.  However,  when  there  is  a 
coupling  between  the  states  and  the  control  inputs,  the  robustness  of  the  system 
is  no  longer  guaranteed.  This  chapter  presents  a  quick  review  of  the  pertinent 
aspects  of  the  standard  theory,  points  out  where  and  how  the  cross  coupling  of 
the  states  and  inputs  occurs,  and  the  problems  this  creates. 

A.   STANDARD  LINEAR  QUADRATIC  CONTROL  THEORY 

A  linear  finite-dimentional  time-invariant  dynamic  system  can  be  described 
in  a  compact  form  by  input  and  output  equations 

x  =  Ax  +  Bu,  (2-1) 

and 

y  =  Oc,  (2-2) 

where  x  is  an  nxl  state  vector,  u  is  an  mxl  control  input  vector,  y  is  a  pxJ 
output  vector,  and  A,  B,  C,  and  D  are  system  matrices  of  the  necessary 
dimensions.   Feedback  control  laws  can  be  used  to  either  stabilize  the  system 
or  modify  its  response  to  meet  the  specifications  of  a  particular  application. 

Assuming  that  all  states  of  the  system  can  be  measured,  a  feedback  control 
law  is  computed  as  a  linear  combination  of  the  states 


u  =  -Kx, 


(2-3) 


where  the  feedback  gain  matrix  K  is  determined  from  the  type  of  problem  and 
the  method  of  control  desired.  This  closed-loop  system  is  depicted  in  Figure  2.1 
Optimal  control  is  a  design  method  which  computes  K  by  minimizing  a  cost 
function  (or  maximizing  a  performance  function).  The  cost  function  is  a  real- 


Figure  2.1   Closed-Loop   Linear  System 

valued  function  of  the  state  of  the  system  and  the  control  input,  both  of  which 
can  vary  with  time.  Cost  functions  take  the  general  form 


J  =  \j{x{t\u{t\t)dt. 


(2-4) 


Application  of  optimal  control  requires  the  specifications  of  the  control  problem 
of  interest  to  be  formulated  in  terms  of  an  appropriate  cost  function.  Minimum 
time  control,  terminal  control,  minimum  control  effort,  minimum  squared  error, 
and  combinations  of  these  are  just  a  few  of  the  possibilities. 

One  cost  function  which  has  been  found  to  be  particularly  useful  is 
quadratic  control  which  takes  the  form 


J  =  j~[x*(t)Qx(t)  +  u(t)Ru(t)\lt,  (2-5) 

where  *  denotes  the  complex  conjugate  transpose,  Q  is  a  positive  semi-definite 
nxn  matrix,  and  R  is  a  positive  definite  mxm  matrix.  Kalman  [Refs.  1&2]  found 
that  the  optimal  feedback  control  law  for  minimizing  this  cost  function  is 

uopt(t)  =  -R-1B*Pxopt(t) 

=  -Kxopt{t\  (2-6) 

where  K  =  R-]B*P  is  determined  by  finding  a  positive  definite  solution  P  for  the 
Steady-State  Riccati  Equation  (SSRE), 

PA  +  JCP  -  PBR~]B*P  +Q  =  0.  (2-7) 

Many  algorithms  are  available  for  finding  P  and  the  feedback  gain  vector  K.  If 
the  system  is  completely  controllable,  P  is  unique  and  the  closed-loop  system 

x  =  {A-BK)x  (2-8) 

will  be  stable. 

B.    ROBUSTNESS  PROPERTIES  OF  LQ  DESIGNS 

To  design  a  control  system,  a  mathematical  model  of  the  system  of  interest 
must  be  developed  so  that  analytical  techniques  may  be  applied.  Linear 
differential  equations  (2-1)  and  (2-2)  comprise  one  such  model.  In  addition,  the 
model  must  be  of  a  form  that  is  accurate  enough  to  give  a  reasonable 
representation  of  the  system,  yet  simple  enough  that  it  leads  to  practical 
analytical  solution  techniques.  Unfortunately,  no  useful  mathematical  model  is 
capable  of  exactly  representing  a  system.   Usually  differences  will  arise  from 


non-linearities  in  a  system  which  is  represented  by  a  linear  model  or  from  the 
use  of  first  or  second  order  models  as  approximations  of  systems  with  higher 
order  dynamics.  In  addition,  parameter  variations  from  operating  conditions, 
environmental  conditions,  measurement  errors,  etc.  may  all  contribute  to 
uncertainties  in  the  operation  of  the  actual  plant. 

It  is  important  to  design  a  control  system  that  can  perform  well  over  as 
much  of  the  expected  range  of  uncertainty  as  possible.  A  system  that  is 
insensitive  to  or  tolerant  of  parameter  variations  and  model  uncertainties  is 
termed  robust.  For  (SISO)  systems,  robustness  is  usually  measured  in  terms  of 
three  parameters  (as  depicted  in  Figure  2.2): 

•  Gain  Margin  —  A  measure  of  how  much  the  open-loop  gain  may  be 

increased  before  the  system  becomes  unstable. 

•  Gain  Reduction  Tolerance  (GRT)  —  A  measure  of  how  much  the  open- 

loop  gain  may  be  reduced  before  the  system  becomes 
unstable. 

•  Phase  margin  ~  A  measure  of  how  much  the  open-loop  phase  delay  may 

be  increased  before  the  system  becomes  unstable. 


Gain  Margin     — 

GRT               1-y 

Phase  Margin 
v^^ NZ           -1 

Irai 

\ 

N^e      \^^^ 

J             *Re 

Figure  2.2  Definitions  of  Robustness  Measures  [After  Ref.  3] 


The  optimal  control  system  shown  in  Figure  2.1   has  an  open-loop 
frequency-domain  transfer  function 

G(s)  =  (sI-AT1B,        s  =  jco  (2-9) 

and  a  loop  gain 


L(s)  =  KG(s) 
:(sl- 
=  R^B^PisI-A^B.  (2-10) 


=  K(sI-Ay1B 


The  robustness  of  this  LQ  system  can  be  determined  by  algebraic  manipulation 
of  (2-7).  First,  add  and  subtract  sP  and  then  regroup  terms: 

sP-sP  +  PA  +  A*P  +  Q-PBR~1B*P  =  0  (2-11) 

-P(sI-A)-(-sI-A*)P  +  Q-PBR~1B*P  =  0.  (2-12) 

Then  multiply  by  G*(s)  =  B*(-sI-A*)-]  from  the  left  and  G(s)  =  (sI-A)-*B  from 
the  right  to  get 

-B\-sl-AY1PB-B*P(sI-Ay1B  +  B*{-sI-A*T1Q{sI-AT1B 

-B*(-sI-A*y1PBR-1B*p(sI-AT1B  =  0  (2-13) 

Noting  that  B*P(sI-A)-1B  =  RL(s),  (2-13)  can  be  rewritten 

-L*{s)R  -  RL(s)  +  G*(s)QG(s)  -  L*  (s)RL(s)  =  0  (2-14) 


or 


[I  +  L(s)]*R[I  +  L(s)]  =  R  +  G*(s)QG(s).  (2-15) 


Since  Q  is  positive  semi-definite  and  G*(s)QG(s)  is  quadratic  in  G(s),  the 
right  hand  side  of  (2-15)  is  greater  than  or  equal  to  R.  Therefore,  the  Nyquist 
plot  of  this  system  must  remain  outside  of  a  unit  circle  centered  at  (-1,  jO)  as 
shown  in  Figure  2.3.  The  optimal  linear  quadratic  controller  is  guaranteed  to 
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Figure  2.3   Robustness  of  an  LQ  System     [After  Ref.  3] 

have  an  infinite  gain  margin,  at  least  60°  of  phase  margin,  and  a  minimum  of 
50%  GRT. 

C.   THE  GENERALIZED  QUADRATIC  CRITERION 

A  more  general  formulation  of  the  quadratic  control  problem  includes 
frequency  shaping  of  the  cost  function.  This  allows  frequency  dependent 
performance  requirements  to  be  included  in  the  system  specifications.  System 
performance  can  be  improved,  but  robustness  may  be  reduced. 


A  frequency  shaped  cost  function  can  be  generated  by  first  transforming 
(2-5)  into  the  frequency  domain  using  Parseval's  theorem, 

J  =  jj™Jx*(jco)QxUcQ)  +  u(jCO)Ru{jCQ)^lCQ.  (2-16) 

The  frequency  shaping  is  accomplished  by  making  Q  and  R  functions  of 
frequency, 

J  =  j  j~o[x*(jQ))QUo))x(jG))  +  u{jco)RUco)u{jco)^ico.  (2-17) 

When  Q(jco)  and  R(jco)  meet  the  additional  constraints  that  they  are  rational 
functions  of  jco  and  decomposable  into 

QUco)  =  fij  +  U\jcD)Q*2  +  Q2U{jco)  +  U*Uco)Q3U(jco),  (2"18) 

where  Yl(jco)  is  a  qxn  matrix  of  rational  functions  of  jco  and  Qj,  Q2,  and  Qs  are 
nxn,  nxq,  and  qxq,  constant  matrices  respectively,  an  optimal  controller  can  be 
obtained  for  this  cost  function  by  augmenting  the  plant  to  eliminate  the 
frequency  dependance  of  the  weighting  matrices  and  then  applying  standard 
LQ  procedures.  [Ref.  4] 

The  augmentation  is  accomplished  by  first  defining 

xUco)  =  n(jco)x(jco)  (2-19) 

and 

u(jco)  =  n(jco)u(jco).  (2-20) 

Then 


ij  -  Dzj  +  Ex, 


(2-21) 


x  —  Lzj  +  Mx, 


(2-22) 


and 


z2  =  FZ2  +  Hu, 


(2-23) 


u  =  Nz2  +  Ju. 


(2-24) 


The  first  term  in  (2-17)  can  now  be  rewritten  using  (2-18),  (2-19),  (2-21),  and 
(2-22)  as 

x*Uco)Q(jco)x(jG))  =  x*(jo))Q1(Jo))xU(o)  +  x*(J(o)Tl*Ua))Qlx(Ja)) 
+x*  (j(D)Q2n{jco)x(jcD)  +  x*  (jo))U*  (jo))Q3n(Jco)x(jco) 


(x\j(D)x\j(D)) 


Qi    Qi 


JV 


rxuay 

*0*o>) 


=  {x  (jco)  ZjiJCO)) 


I    M 

o    £ 


Qi    Qi 
Q*2   Qs 


i    o' 

M    L 


rx{jcoY 


(2-25) 


Similarly,  the  second  term  becomes: 


u  {j(D)R{jco)u(jQ))  =  {z2{jco)  u  (jco)) 


N*    0 
/     I 


Rj     R2 

lR2      R3j 


N    J 
0     I 


rZ2(jCD)' 

\u(jco) 


(2-26) 


Finally,  the  whole  system  and  the  cost  function  can  be  written  as: 
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f  +  \ 


A     0     0 
EDO 


{z2J    10     0     F]{z2) 


fx^ 


Zl 


0 


u 


=  A 


'^ 


K.Z2) 


+  Bu 


(2-27) 


and 


J  =  i\_oo(x  *  z)  £  M*) 


rj    r2 
r:   n 


o 

0 


0 
0 


0 

0' 

(X) 

0 

0 

Zl 

*i 

% 

*2 

r2 

Vs. 

KUJ 

dco 


i  r°°  /    *  *    *    *x 
=  l]_00(x    Zj  z2  u  ) 


Q    s 

S*     R 


(x  zj  z2  u)  dco, 


(2-28) 


where 


1}  =  Qj  +  M  Q2  +  Q2M  +  M  Q3M, 

r2  =  Q2L  +  M*Q3L, 

r3  =  L*Q3L, 

f7  =  N*R3N, 

tF2  =  N*R3J  +  N*R2, 

XF3=R1  +  R2J  +  J*R*2  +  J*R3J, 


Tj 

r2    o' 

Q  = 

r; 

r3    o 

0 

0     Yj 

s*  =  (o 

o  r2), 

R  = 

*S- 
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Note  this  state  augmentation  procedure  has  generated  a  cross-state-input 
matrix  S  which  has  serious  implications  for  the  robustness  of  this  control 
system.  The  gain  and  phase  margins  and  gain  reduction  tolerance  of  the 
standard  LQ  system  are  no  longer  guaranteed. 

D.   ROBUSTNESS  FOR  A  GENERALIZED  QUADRATIC 
CRITERION 

The  addition  of  frequency  shaping  in  the  cost  function  generated  a  cross - 

state-input  term.   A  procedure  very  similar  to  those  performed  on  (2-9)  through 

(2-15)  is  developed  below  to  identify  the  effect  of  S  on  the  robustness  of  the 

system.  First,  the  following  definitions  are  made  to  simplify  the  expression: 


JC  = 


x 

lZ2] 


(2-29) 


Furthermore,  define 


-i , 


u,  =  -R'1B  Px, 


(2-30) 


and 


u  =  u'  -R  *S*x 


(2-31) 


Then 


x  =  Ax  +  Bu 

=  Ax  +  B(u'-R~1S*x) 

=  {A-BR-1S*)x^Buf 
=  A'x  +  Bu'. 


(2-32) 
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From  (2-28), 

*  *  *      *  * 

x  Qx  +  x  Su  +  u  S  x  +  u  Ru 

=  xQx  +  x*S(u'  -  R'JS*x) +(u'-  RT1S*x)*  S*x 

+("'  -  R~JS*x)*  R{u'  -  R^S'x) 
=  x*(Q-SR~1S*)x  +  u'*Ru 
=  x*Q'x  +  u'*Ru.  (2-33) 

Using  the  above  definitions,  the  SSRE  becomes 

PA'  +  A'*P  +  Q'-  PBR~1B*P 

=  P{A-BR-1S*)  +  (A-BRT1S*)*P+Q 


-SR-'S*  -PBR^B'P 


*-i 


=  PA  +  A  P+Q-(PB+S)R~1(PB+S)  =0.  (2-34) 


Finally, 


7      * 

u  =  u'  -R   S  x 
=  -R~1B*Px-R'1S*x 
=  -R~1(B*P+S*)x,  (2-35) 


and 


K  =  R-1(B*P+S*)  =  R-\PB  +  S)*  (2-36) 

Now,  given  the  modifications  as  in  (2-29)  through  (2-36)  the  analysis  of  the 
system  proceeds  just  as  in  the  standard  LQ  case: 

G(s)  =  (sI-AT1B,  (2-37) 
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and 


L(s)  =  KG(s)  =  K(sI-AT1B.  (2-38) 


Using  (2-36), 


L(s)  =  R~1(PB  +  S)\sI-AT1B.  (2-39) 

Adding  and  subtracting  sP  to  the  SSRE  and  rearranging  yields 

-P(sI-A)-(sI-A*)P+Q-(PB+S)R-1(PB  +  S)*  =  0.  (2-40) 

Multiplying  from  the  right  by  G(s)  and  from  the  left  by  G  *(s)  and  rearranging 
terms  as  was  done  in  (2-14)  and  (2-15)  results  in 

[I  +  L(s)]*R[I+L(s)]  =  R+G*(s)QG(s)  +  S*G(s)  +  G*(s)S.        (2-41) 

Comparing  (2-41)  with  (2-15),  there  are  two  new  terms  as  a  result  of  the 
cross- state-input  matrix  S.  The  standard  optimal  LQ  system  was  guaranteed 
an  infinite  gain  margin,  50%  GRT,  and  60°  of  phase  margin.  If  these  two  new 
terms  in  (2-41)  result  in  reducing  the  magnitude  of  the  right  hand  side,  they  will 
effectively  reduce  the  robustness  of  the  system.  Of  course,  the  converse  may 
also  be  true.  In  some  systems,  the  right  hand  side  may  be  increased  and 
robustness  will  correspondingly  increase.  The  effects  of  S  will  depend  on  the 
particular  system  and  the  required  frequency  shaping  in  the  cost  function. 
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III.     MODEL   UNCERTAINTY   AND   ROBUSTNESS 
MEASURES  FOR  MULTIVARIABLE  SYSTEMS 

Multi-input  multi-output  (MIMO)  systems  introduce  significantly  greater 
complexity  into  the  control  design  problem.  The  additional  inputs  and  outputs 
provide  a  designer  with  much  more  flexibility  in  the  solution  to  the  problem  of 
stabilizing  a  system,  but  they  also  complicate  the  subject  of  defining  the 
robustness  of  the  system.  This  chapter  presents  some  concepts  used  for 
MEMO  systems  and  develops  the  robustness  measure  to  be  used  for  this  thesis. 

A.  MODEL  UNCERTAINTY  FOR  MIMO  SYSTEMS 

Since  no  design  model  will  exactly  match  the  physical  plant,  it  is  important 
to  consider  these  model  uncertainties  for  control  system  design.  The  form  used 
in  representing  them  depends  on  the  degree  to  which  the  physical  mechanisms 
of  the  uncertainties  are  understood  or  known.  In  general,  representations  can 
be  placed  into  two  categories  based  on  the  amount  of  structure  in  the 
uncertainty.  Highly-structured  representations  are  those  in  which  the  amount  of 
parameter  variation  can  be  confined  within  well  defined  limits  for  models  at 
various  operating  points.  An  example  would  be  the  aerodynamic  coefficients  of 
a  missile  for  different  altitudes  and  velocities.  Less -structured  representations 
usually  place  bounds  on  the  uncertainty  of  the  plant  transfer  function.  Doyle 
[Ref.  5]  gives  two  general  ways  of  representing  this: 

G\jco)  =  GUco)  +  AGUco),  (3-1) 

where 
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o[AG(jco)]<la(CQ)  forVco>0  (3-2) 


and 


G\j(0)  =  G{jco)[l  +  A(jco)]t  (3-3) 


where 


o[A{j(0)]<lm{(0)  forVco>0.  (3-4) 

The  maximum  singular  value  a  is  defined  by 

a[A]  =  max||Ax||  =  JkJ^a]  (3-5) 

Neither  (3-1)  and  (3-2)  nor  (3-3)  and  (3-4)  attempt  to  give  any  indication  of  the 
causes  of  the  uncertainties.  The  additive  uncertainty  representation,  (3-1)  and 
(3-2),  bounds  the  matrix  G'  within  a  region  la.  about  G  The  multiplicative 
uncertainty  representation,  (3.3)  and  (3-4)  (which  is  preferred  here  since  it 
applies  the  same  uncertainty  to  both  uncompensated  and  compensated  transfer 
functions),  confines  the  matrix  G'  to  a  normalized  neighborhood,  /m>  about  G 
This  type  of  uncertainty  can  be  incorporated  into  the  system  model  as  shown  in 
Figure  3.1.  Knowledge  about  and  experience  with  the  uncertainties  affecting  a 
particular  system  are  the  determining  factors  in  the  magnitude  of  lm. 
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Figure  3.1   System  with  Unstructured  Uncertainties 

B.  ROBUSTNESS  MEASURES  FOR  MIMO  SYSTEMS 

Many  methods  have  been  investigated  and  developed  for  measuring  the 
robustness  of  MIMO  systems.  The  principal  region  and  minimum  singular 
value  techniques  will  be  discussed  here. 

1.     Principal  Region 

Postlethwaite  [Ref  6]  has  developed  a  method  of  extending  the  gain 
and  phase  margin  concept  to  MIMO  systems  through  the  definition  of  a 
principal  region.  The  principal  region  is  found  by  polar  decomposition  of  the 
loop  transfer  matrix  L(s): 


L(s)  =  U(s)HR(s) 


(3-6) 


and 


L(s)  =  HL(s)U(s), 


(3-7) 


where  U(s)  is  unitary  and  HR(s)  and  HL(s)  are  positive  semi-definite  Hermitian 
matrices  with  the  same  eigenvalues.    Principal  gains  (or  singular  values)  are 
defined  as  the  eigenvalues  of  HR(s)  or  HL(s)    and  the  principal  phases  are 
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defined  as  the  phases  of  the  eigenvalues  of  U(s).  The  maximum  and  minimum 
principal  gains,  omax  and  om[n,  and  the  maximum  and  minimum  principal 
phases,  fymax  and  tymin-,  for  each  s,  form  a  curvilinear  rectangle  in  the  Nyquist 
plane.  Combining  all  the  rectangles  for  values  of  s  on  the  Nyquist  D  contour 
forms  the  principal  region.  A  gain  margin,  phase  margin,  and  gain  reduction 
tolerance  can  be  defined,  as  depicted  in  Figure  3.2,  at  the  boundary  of  the 
principal  region. 


Im 


Re 


Principal  Region 


Figure  3.2  Principal  Region  and  Stability  Margins  for  MIMO 

systems  [After  Ref.  4] 
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This  method  generates  graphical  information  which  does  not  easily 
lead  to  defining  a  single  parameter  which  may  be  used  to  adjust  performance 
and  robustness.  However,  the  use  of  just  the  minimum  singular  values,  while  a 
more  conservative  approach,  provides  similar  robustness  information  and  does 
lend  itself  to  the  goal  of  this  thesis. 
2.     Minimum  Singular  Value 

A  measure  of  the  robustness  of  a  system  with  the  uncertainties 
described  above  can  be  obtained  as  follows: 


and 


Then 


and 


z  =  KG(I  +  A)e  (3-8) 


e  =  r-z.  (3-9) 


e  =  r-KG{I  +  A)e,  (3-10) 


(I  +  KG(I  +  A))e  =  r,  (3-11) 


e  =  (I  +  KG{I  +  A)T1r 
=  (I  +  KG  +  KGA)~1r 

=  [{I  +  KG){I  +  {I  +  KG)'1  KGA)]'1  r 

=  [l  +  {I  +  KGT1KGA]~\l  +  KG)~1r.  (3-12) 
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Using  the  small  gain  theorem  and  (3-12),  a  sufficient  condition  to 
guarantee  closed-loop  stability  is 


m  (I  +  KGy'KG  \[\\  A    <1  (3-13) 

or  equivalently 

a[  (I  +  KG)'1  KG  ]a[  A  ]<i.  (3-14) 

Using  (3-4), 

a[  (I  +  KG^KG  ]<  —  (3-15) 

''tn 

Equation  (3-15)  can  be  written  in  a  more  convenient  form  using  some  results  of 
matrix  theory.  First, 

(I  +  KG^KG^V  +  iKGT1)-1,  (3-16) 

so  (3-15)  becomes 

a[(/+(*G)-V]</-.  (3-17) 


and,  using 


o[A]' 
equation  (3-17)  becomes 


oM  =  -T77.  (3-18) 
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all  +  iKGT1]     lm 


or 


o[l  +  (KGy1]>lm.  (3-20) 

Therefore,  the  robustness  measure  to  be  used  here  is  defined  as 

J^all  +  iKGT1],  (3-21) 

i.e.,  Jr  is  the  minimum  singular  value  for  the  system  in  the  frequency  range  of 
interest.  Jr  is  a  single  number  which  gives  a  measure  of  the  robustness  of  the 
system.  For  the  standard  optimal  LQ  system,  Jr  is  guaranteed  to  be  >  0.5. 
Therefore,  the  robustness  of  a  system  can  be  computed  and  directly  compared 
to  the  standard  optimal  system  and  a  determination  made  as  whether  or  not  it  is 
satisfactory. 

As  stated  before,  this  method  is  very  conservative.  There  are  the 
same  number  of  singular  values  as  the  dimension  of  the  system  and  it  is  unlikely 
that  the  uncertainties  in  the  model  will  affect  the  whole  system  equally  across 
all  the  states.  Therefore,  some  models  may  very  well  have  much  more 
robustness  than  what  is  predicted  by  this  method. 
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IV.     ROBUST   SUBOPTIMAL   DESIGN 

The  general  idea  for  the  suboptimal  robust  design  is  to  trade  off  optimality 
(i.e.,  performance)  for  a  greater  degree  of  robustness.  In  the  case  of  the 
generalized  linear  quadratic  control  system,  this  means  computing  the  feedback 
gains,  K  =  R~\PB  +  S)* ,  by  modifying  Q,  R,  or  S  in  (2-41). 

A  number  of  design  possibilities  exist  and  seven  options  investigated  under 
this  thesis  work  are  presented  below. 

A.  ROBUSTNESS   DESIGN   OPTIONS 

Equation  (2-41)  is  repeated  here  as  it  serves  as  the  basis  for  evaluating 
robustness: 

(I  +  L)*R(I  +  L)  =  R  +  G*QG  +  S*G  +  G*S,  (4-1) 

where  L  and  G  are  understood  to  be  functions  of  s  =jco.   In  addition,  using  the 
relation 

(R  +  S*G)*  R~'(R  +  S*G)  =  R  +  S*G  +  G*S  +  G*SR-*S*G,  (4-2) 

equation  (4-1)  can  be  rewritten  as 

(I  +  L)*R(I  +  L)  =  G*(Q-SR-1S*)G  +  (R  +  S*G)*R-1(R  +  S*G).      (4-3) 

Note  that  the  right  side  is  guaranteed  to  be  greater  than  or  equal  to  zero  if 

Q-SR-JS*>0.  (4-4) 
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The  greater  the  right  side  of  (4-3)  the  better  the  robustness  and  it  certainly 
needs  to  be  greater  than  zero  to  guarantee  stability.  All  options  are  presented 
using  a  tilde,  ~,  to  indicate  the  modified  parameters. 

1.  Option  I:    S  =  0 

The  first  option  is  to  simply  ignore  the  cross-state-input  matrix  by 
setting  S  =  S  =  0.  Equation  (4-1)  would  be  reduced  to 

(I  +  L)*R(I  +  L)  =  R  +  G*QG,  (4_5) 

which  is  the  same  as  (2-15)  for  the  standard  LQ  situation.  However,  this  is  just 
a  special  case  of  option  IV. 

2.  Option  II:    Q=Q  +  SS\  R  =  R  +  I 
These  substitutions  yield 

(/  +  L)*  R(I  +  L)  =  R  +  G*QG  +  S*G  +  G*S  (4-6) 

or 

(I  +  L)*(R  +  I)(I  +  L)  =  R  +  G*QG  +  (I  +  S*G)\l  +  S*G).  (4-7) 

Here  the  terms  on  the  right  hand  side  have  either  been  forced  into  quadratic 
forms  or  are  positive,  ensuring  that  robustness  is  improved. 

3.  Option  III:    R  =  pR 

This  is  a  straight-forward  introduction  of  a  multiplicative  parameter 
into  (4-1)  which  gives 

(I  +  L)*R(I  +  L)  =  R  +  G*QG  +  S*G  +  G*S 

=  pR  +  G*QG  +  S*G  +  G*S  (4-8) 
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or 

(I+L)*R(I  +  L)  =  R  +  -(G*QG  +  S*G  +  G*S).  (4-9) 

P 

As  p  is  increased,  the  right  side  approaches  R.  If  the  robustness  of  the 
system  had  originally  been  degraded  by  S,  theoretically,  the  best  that  could  be 
achieved  by  this  option  is  exactly  the  minimum  guaranteed  by  the  standard  LQ 
system. 

4.  Option  IV:    S  =  pS,   p<l 

As  in  option  HI,  this  is  also  a  simple  substitution  using  a  multiplicative 
parameter.  This  time,  however,  it  is  applied  direcdy  to  S: 

(I  +  L)*R(I  +  L)  =  R+G*QG  +  S*G  +  G*S 

=  R  +  G*QG  +  p(S*G  +  G*S)  (4-10) 

Note  that  when  p  =  0,  (4-9)  reduces  to 

(I  +  L)*R(I  +  L)  =  R+G*QG,  (4-11) 

the  standard  optimal  LQ  case  and  option  I  is  therefore,  a  special  subset  of  this 
option. 

5.  Option  V:    R  =  pR,  Q  =  Q  +  (p-l)SR-1S\  p>l 

Substituting  into  (4-3): 

(I  +  L)*R(I  +  L)  =  G*[Q-SR-1S*]G  +  (R  +  S*G)*R-1(R+S*G).  (4-12) 
This  results  in 


24 


(/  +  L)*  pR(I  +  L)  =  G*[Q  +  (p-  2)SR~1S*]G 

+(pR  +  S*G)\pRy1(pR  +  S*G)  (4-13) 


or 


(I  +  L)*R(I  +  L)  =  -G*[Q  +  (p-2)SR-1S*]G 
P 

+(R  +  -S*G)*  R~1(R  +  -S*G).  (4-14) 

P  P 

As  p  is  increased,  (4-14)  approaches 

(I+L)*R(I+L)  =  R  +  G*SR~1S*G.  (4-15) 

Again,  the  term  with  S  has  been  made  quadratic  in  order  to  improve 
robustness. 

6.     Option  VI:    Q  =  pQ 

Here,  a  simple  substitution  of  a  multiplicative  parameter  into  (4-3): 

(I  +  L)*R(I+L)  =  G*[Q-SR~1S*]G 

+(R  +  S*G)*R-1(R+S*G)  (4-16) 

or 

rV; 

HR  +  S'GfR-'iR  +  S'G).  (4- 


(I  +  L)  R(I  +  L)  =  G  [pQ-SRT'S  ]G 


Note  that  the  right  side  is  guaranteed  to  be  greater  than  or  equal  to  zero  only  if 

pQ-SR-]S*>0.  (4-18) 
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For  optimal  LQ  theory,  Q  is  only  required  to  be  positive  semi-definite.  For  this 
option,  there  is  no  guarantee  that  pQ-SR^S*  can  be  made  greater  than  zero 
unless  Q  is  positive  definite. 

7.     Option  VII:    Q  =  Q  +  pSR_1S%  p>l 
Substituting  Q  into  (4-3)  yields: 

(I+L)*R(I  +  L)  =  G*[Q-SR~1S*]G 

+(R  +  S*G)'¥R~1(R+S*G) 

=  G*[Q  +  (p-l)SR~1S*]G 
HR+S*G)*R-](R+S*G)  (4-19) 

Again,  the  right  side  is  greater  than  or  equal  to  zero  for  p  >  1. 

This  list  certainly  does  not  exhaust  the  possibilities,  but  it  does  provide  many 
avenues  for  investigation  of  the  robustness  of  an  MIMO  system  optimized  by  a 
general  quadratic  criterion.  It  should  also  be  noted  that  when  Q,  R,  or  S  is 
modified,  the  control  system  is  no  longer  optimal. 
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V.  A   NUMERICAL   EXAMPLE   AND   SIMULATION 

To  explore  the  utility  of  the  theory  so  far  developed,  a  numerical  example 
will  be  presented  and  the  control  systems  using  the  various  options  in  the 
preceding  chapter  will  be  designed  and  simulated.  The  model  to  be  used  is  the 
one  used  by  Chen  [Ref.  3].  The  development  and  explanation  of  this  model  is 
briefly  summarized  here.  The  simulations  will  use  a  wind  gust  as  the 
disturbance  input  to  the  system  to  test  the  transient  response  of  various 
suboptimal  designs. 

A.  STATICALLY  UNSTABLE  AIRCRAFT  MODEL 

The  model  for  this  simulation  is  the  F-20  Tigershark  designed  by  Northrop 
Corp.  The  aircraft  is  designed  to  be  extremely  agile,  and  as  a  result,  is  highly 
unstable  in  the  absence  of  a  control  system.  In  addition,  the  simulation  will  be 
based  on  a  special  flight  mode  called  fuselage-pitch-pointing  (FPP)  mode.  The 
fuselage-pitch-pointing  concept  provides  the  ability  to  alter  the  pitch  without 
changing  the  flight  path  (Figure  5.1).  Another  interesting  flight  mode  called 
direct  lifting  (DL)  is  shown  in  Figure  5.2,  where  the  aircraft's  altitude  can  be 
changed  without  changing  the  pitch. 

These  are  not  the  intuitive  ways  of  flying  an  aircraft  and  therefore,  would 
not  be  the  normal  modes  of  operation.  Instead,  it  might  be  used  in  special 
circumstances  such  as  bombing  runs  or  during  dogfights  where  the  ability  to 
rapidly  change  altitude  or  pitch  of  the  aircraft  independently  may  provide 
significant  advantage  over  an  adversary.  In  addition,  the  FPP  mode  provides  a 
good  MEMO  model  for  investigating  robust,  suboptimal  controller  design. 
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Right  Path 


Figure  5.1   Pitch  Pointing,  Constant  Flight  Path 


Figure  5.2   Direct  Lifting,  Constant  Pitch 

A  linearized  pitch-axis  model  for  the  F-20  at  Mach  0.5  and  an  altitude  of 
10,000  feet  is: 
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The  state  variables  are  angle  of  attack  a,  pitch  rate  q  and  pitch  angle  6,  as 
shown  in  Figure  5.3. 


Figure  5.3   Definitions  of  State  Variables  a,  q,  and  8 

The  two  control  inputs  are  stabilizer  position  8h  and  trailing  edge  flap  position 
5f.  This  model  has  three  poles  located  at  s  =  0,  -2.519,  and  0.795. 

The  design  problem  is  to  develop  an  FPP  control  system  that  will  reject 
turbulence  and  yet  provide  good  robustness.  This  thesis  will  focus  on  a 
regulator  type  of  control  which  will  endeavor  to  maintain  the  states  of  the 
system  at  desired  values  in  the  presence  of  disturbances  to  the  system.  An 
appropriate  frequency-shaped  quadratic  cost  function  which  accomplishes 
these  objectives  is 


Jc  =j~Jl  50a2 +  q2  + 50(6 -a)  + 
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Turbulence  is  typically  manifested  by  disturbance  at  frequencies  <  30  rad/sec. 
Therefore,  the  cost  function  is  designed  with  a  high  pass  filter  on  the  control 
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inputs  to  allow  greater  control  effort  at  low  frequencies  and  to  penalize  control 
effort  at  higher  frequencies.  In  addition,  flight  path  angle  y  =  6  -  a  is  heavily 
weighted  to  limit  deviations  from  the  flight  path.  Turbulance  however,  is  a 
transient  phenomenon,  therefore  this  optimization  focuses  more  on  rejecting 
these  transitory  disturbances  in  pitch  rate,  pitch,  and  angle  of  attack  than  on 
steady-state  deviations  from  the  flight  path. 

For  this  problem,  the  frequency  dependent  weighting  matrix  on  the  control 
is  selected  to  be 
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and  it  can  be  readily  decomposed  into 
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which  meets  the  criteria  of  (2-22). 

Following  the  augmentation  procedure  described  in  Chapter  II,  this  system 
can  be  expressed  in  a  frequency-independent  format.  First,  expand  the  input 
terms  of  (5-3)  into  the  form: 
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Now  define 
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and  let 
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This  results  in  the  new  state  equations: 


z1  =  -30z]  +  30Sh 


(5-9) 


and 
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(5-10) 


Similar  steps  for  8h  result  in: 


z2  =  -30z2  +  30Sf 


(5-11) 


and 


u2  =  -z2  +  8f, 


(5-12) 


Thus  the  new  augmented  system  is 
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and  the  frequency-independent  cost  function  is 
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B.   SIMULATIONS 

The  simulations  will  use  a  vertical  wind  gust  with  a  maximum  velocity  of 
62.5  feet  per  second.  This  disturbance  S  will  be  introduced  into  the  system  as 
shown  in  Figure  5.4.  The  wind  gust  effectively  changes  the  direction  of  the 
incident  wind  and  consequently  the  angle  of  attack  of  the  aircraft,  as  indicated 
in  Figure  5.5. 


Figure  5.4   Wind  Gust  Disturbance  Inputs 


Flight  path 


Wind  gust 


Figure  5.5  Angle  of  Attack 
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All  performance  Jc  and  robustness  Jr  values  for  theses  simulations  are 
computed  using  the  rectangular  wind  gust.  For  comparison  purposes,  the 
optimal  system  response  is  presented  first. 

1.     Optimal  System 

Optimization  of  this  control  problem  results  in 

K  = 


-11.1747    -1.2199    -0.5267    -0.3593    -0.0169 
-37.0471     3.1484     35.2571    -0.6706    -0.5002 


Jc  =  3.9678, 

Jr  =  0.3529, 
and  Figure  5.6  shows  the  minimum  singular  values,  as  defined  in  (3-21),  as  a 
function  of  frequency.  Note  that  within  the  range  of  frequencies  expected  for 
turbulence  (<  30  rad/sec),  the  plot  drops  below  0.5.  Figures  5.7  and  5.8  show 
that  the  response  to  turbulence  is  quite  good.  However,  with  a  robustness  of 
0.3529,  GRT  is  only  35  percent,  leaving  a  much  reduced  margin  for  errors  or 
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Figure  5.6   Minimum  Singular  Values  of  the  Optimal  System 
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Figure  5.7  Response  of  the  Optimal  System  to  a  Rectangular  Gust 
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uncertainties.  A  suboptimal  design  will  be  able  to  trade  off  some  of  the  desired 
performance  in  order  to  recover  greater  tolerance  of  model  and  parameter 
error. 

2.     System  with  Q  =  Q  +  SS*,  R  =  R  +  I 

This  simple  modification  to  Q  and  R  results  in  the  singular  values 
shown  in  Figure  5.9,  with 
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Jc=  11.1382, 
Jr  =  0.8049. 


1 

6 

5 
4 
3 
2 
1- 

0 

] 

j...;...i.i.j.j.;;j 

11 

■■:-*-!■: 

LO"2                        10-1                        10°                         101                          102 

co  (rad/sec) 

Figure  5.9   Minimum  Singular  Values  Using  Q  =  Q  +  SS%  R  =  R  +  I 

This  is  a  significant  increase  in  robustness.  Unfortunately,  it  also  produces  a 
significant  loss  of  performance,  as  depicted  in  Figures  5.10  and  5.11.  There  is 
nearly  an  order  of  magnitude  increase  in  pitch  angle,  at  least  twice  as  much 
pitch  rate,  and  a  significant  drop  in  response  time.    Additionally,  there  is  no 
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adjustable  parameter  available  for  selecting  some  intermediate  point  between 
this  and  the  optimal  system. 
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Figure  5.10    Square  Gust  Response  of  Q  =  Q  +  SS  ,  R  =  R  +  I 
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3.     System  with  R  =  pR 

This  option  provides  for  adjustment  of  the  amount  of  trade  off.  As  can 
be  seen  in  Figures  5.12  and  5.13,  some  improvement  in  robustness  is  obtained, 
but  the  maximum  for  this  system  is  only  0.4583  for  p  =  3.4.  This  option  appears 
to  provide  very  limited  and  insufficient  improvement. 
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Figure  5.12    Performance  vs  p  for  R  =  pR 


Figure  5.13    Robustness  vs  p  for  R  =  pR 


38 


4.     System  with  S  =  pS,  p<l 

This  option  appears  to  hold  great  promise.  Comparing  Figures  5.14 
and  5.15  with  5.12  and  5.13,  this  scheme  provides  significantly  greater  increase 
in  robustness  for  a  much  smaller  loss  of  performance.  Moreover,  the  use  of  a 
simple  multiplicative  parameter  p  allows  control  engineers  to  adjust  the  amount 
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engineers  to  adjust  the  amount  of  trade-off.   Figures  5.16  and  5.17  demonstrate 
the  response  to  the  two  disturbances. 
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Figure  5.17    Triangular  Gust  Response  for  S  =  pS,   p<l 
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5.     System  with  R  =  pR,  Q  =  Q  +  (p-l)SR_1S*,  p>l 

Figures  5.18  and  5.19  show  that  although  Jr  does  get  above  0.5,  there 
is  still  limited  improvement  for  the  increased  cost.  Apparently  applying  the 
modifications  directly  to  the  control  weighting  matrix  provides  little  benefit.  The 
responses  to  the  disturbances  are  presented  in  Figures  5.20  and  5.21. 
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Figure  5.18    Performance  vs  p  for  R  =  pR,  Q  =  Q  +  (p-l)SR_1S  ,   p>l 
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Figure  5.19    Robustness  vs  p  for  R  =  pR,  Q  =  Q  +  (p-l)SR_1S  ,  p>l 
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6.     System  with  Q  =  pQ 

From  Figures  5.22  and  5.23  it  appears  that  the  cost  is  increasing  nearly 
linearly  with  p.  However,  for  1  <  p  <  5  the  robustness  increases  dramatically 
while  the  cost  only  goes  up  to  about  25.   This  isn't  as  good  as  using  option  IV, 
where  the  cost  increased  far  less  for  a  comparable  increase  in  robustness,  but  it 
does  show  increased  robustness  with  a  variable  parameter.   Figures  5.24  and 


Figure  5.22    Performance  vs  p  for  Q  =  pQ 


Figure  5.23    Robustness  vs  p  for  Q  =  pQ 
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5.25  show  that  the  responses  for  both  types  of  disturbance  are  still  quite  good. 
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Figure  5.24    Square  Gust  Response  for  Q  =  pQ 
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Figure  5.25    Triangular  Gust  Response  for  Q  =  pQ 
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7.     System  with  Q  =  Q  +  pSR_1S*,  p>l 

Figures  5.26  through  5.29  show  that  this  option  gives  excellent  results, 
similar  to  and  slightly  better  than  option  IV.  Using  p  =  10  gives  a  robustness  of 
nearly  0.8  and  increases  the  cost  to  only  about  11.7.  It  would  appear  that 
options  IV  and  VII  would  be  the  best  choices  for  suboptimal  robust  control 
design. 
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Figure  5.26    Performance  vs  p  for  Q  =  Q  +  pSR_1S  ,  p>l 
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Figure  5.27    Robustness  vs  p  for  Q  =  Q  +  pSR  XS  ,  p>l 
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C.   FURTHER  INVESTIGATION  OF  OPTION  IV 

As  shown  in  the  previous  section,  this  suboptimal  design  displayed 
significant  promise.  In  this  section,  the  effects  of  p  are  further  studied. 
Equation  (2-32)  is  repeated  here: 


T  7    f°°     /         *       *         *  *\ 


r7 

r2 

0 

0' 

(X) 

r; 

r3 

0 

0 

Zl 

0 

0 

% 

^2 

z2 

0 

0 

% 

*J. 

<UJ 

~Q 

s' 

R 

(x  Zj 

z2  u) 

Tdco 

dco 


(5-15) 


It  is  evident  that  multiplying  S  by  p  is  in  effect  multiplying  *F2  by  p.  But, 


Y2  =  N  R3J  +  N  R2, 


(5-16) 


and  therefore 


pV2=p(N  R3J  +  N  R2) 
=  pN*R3J  +  pN*R*2 
=  N*R3J  +  N*R2. 


(5-17) 


Since  N  is  produced  by  the  augmentation  of  the  inputs  as  given  in  (2-28),  the 
physical  result  then,  is  a  modification  of  Tl(jco),  the  frequency  shaping  of  the 
control  inputs. 

For  the  aircraft  system  in  this  example,  (5-10)  and  (5-12)  give 


uj  =  Nzj  +  J8h 
=  -zi  +  $h 


(5-18) 
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and 


u2  =  Nz2  +  JSf 
=  -z2  +  Sf. 


(5-19) 


Working  backwards  through  the  augmentation  procedure, 


Uj  = 


30N^ 


s  +  30 


Sh 


(s  +  30)J-30N^ 


s  +  30 


Sh 


s  +  30- 


30(p)\ 


s  +  30       J 


(5-20) 


and 


u2  = 


fs+30 


30(p)\ 


K       s  +  30       J 


'/• 


(5-21) 


For  p  =  l,  (5-20)  and  (5-21)  reduce  to  the  original  high-pass  filter  and  the 
resulting  response  and  robustness  are  those  of  the  optimal  system.  Decreasing 
p  to  -1  in  several  steps  alters  the  frequency  shaping  as  shown  in  Figure  5.30. 
Evidently,  this  aircraft  system  is  more  robust  if  low-pass  filtering  is  used  on  the 
control  inputs. 

Another  interesting  point  is  that  p  had  to  be  limited  to  [-1,  1]  for  this 
system.  Outside  these  limits,  Q^Q-SR^S*  was  no  longer  positive  semi- 
definite  as  required  by  optimal  control  theory  to  ensure  a  stable  close-loop 
system. 
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VL     CONCLUSION 

The  frequency  shaping  in  the  generalized  quadratic  cost  function  can  be 
converted  to  the  standard  linear  quadratic  control  form  by  state  augmentation. 
However,  this  process  generates  coupling  between  the  states  and  the  control 
inputs  and  eliminates  the  guaranteed  good  robustness  characteristics  of  LQ 
optimal  designs.  This  thesis  has  explored  and  developed  a  design  technique  for 
improving  the  robustness  of  these  systems.  This  technique  involves 
modification  of  the  stability  equations  to  regain  robustness  at  the  expense  of 
optimality.  The  following  modifications  were  studied: 

•  Option  I:  S  =  0 

•Optionll:  Q  =  Q  +  SS*,  R  =  R  +  I 

•  Option  IH:  R  =  pR 

•  Option IY:  S  =  pS,  p<l 

•  Option  V:  R  =  pR,  Q  =  Q  +  (p-l)SR_1S\  p>l 

•  Option  VI:  Q  =  pQ 

•  Option  VH:  Q  =  Q  +  pSR_1S* ,  p  >  1 

All  of  these  involve  modifications  to  the  weighting  matrices  in  the  cost  function. 
The  unstable  aircraft  model  used  in  the  simulations  in  this  thesis  responded  best 
to  Options  IV  and  VH. 

This  thesis  provides  control  systems  designers  with  the  ability  to  perform  a 
trade-off  between  optimized  performance  and  robustness.  This  is  accomplished 
through  the  introduction  of  an  adjustable  parameter  in  the  governing  stability 
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equations.    This  approach  provides  a  direct  and  straight-forward  tool  which 
designers  can  easily  apply  to  their  specific  control  problems. 
In  addition,  some  areas  identified  for  future  research  include: 

•  Analysis  of  the  physical  meaning  of  the  changes  made  by  each  of  the 

design  options,  similar  to  that  done  in  Chapter  V  for  Option 
IV. 

•  More  comprehensive  comparison  of  all  the  options  with  each  other  to 

identify  their  similarities  and  differences  in  response  for 
various  values  of  p. 

•  Study  of  the  effectiveness  of  this  technique  for  systems  in  which 

observers  or  Kalman  filters  are  necessary  due  to  the  use  of 
partial  state  feedback. 


55 


REFERENCES 


[1]  R.  E.  Kalman,  "Contributions  to  the  Theory  of  Optimal  Control, "  Pol.  Soc. 
Mat.  Mex.,  Vol.  5,  pp  102-119, 1960. 

[2]  R.  E.  Kalman,  "When  is  a  Linear  System  Optimal?"  Trans.  ASME,  J. 
Basic  Eng.,  Vol.  86,  pp  1-10, 1964. 

[3]  W.  Z.  Chen,  "Constrained  Optimization  Approach  to  Feedback  Control  of 
State  Augmented  Systems, "  Ph.D.  Dissertation,  University  of  California, 
Los  Angeles,  California,  1990. 

[4]  W.  Z.  Chen,  "Parameter  Optimization  for  an  H2  Problem  with 
Multivariable  Gain  and  Phase  Margin  Constraints, "  presented  at  the 
American  Controls  Conference,  Chicago,  Illinois,  24-26  June  1992. 

[5]  J.  C.  Doyle  and  G.  Stein,  "Multivariable  Feedback  Design:  Concepts  for 
a  Classical/Modern  Synthesis"  IEEE  Trans,  on  Automatic  Control,  Vol 
AC-26,No.  l,pp  4-16,  1981. 

[6]  I.  Postlethwaite,  J.  M.  Edmunds,  A.  G.  J.  MacFarlane,  "Principle  Gains 
and  Principle  Phases  in  the  Analysis  of  Linear  Multivariable  Feedback 
Systems, "  IEEE  Trans,  on  Automatic  Control,  Vol.  AC-26,  No.  1,  pp  32- 
46, 1981. 


56 


INITIAL   DISTRIBUTION   LIST 


No.  Copies 


1.  Defense  Technical  Information  Center  2 
Cameron  Station 

Alexandria,  Virginia  22314-6145 

2.  Library,  Code  52  2 
Naval  Postgraduate  School 

Monterey,  California  93943-5000 

3.  Chairman,  Code  EC  1 
Department  of  Electrical  &  Computer  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

4.  Prof.  Won-Zon  Chen,  Code  EC/Cw  1 
Department  of  Electrical  &  Computer  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

5.  Prof.  Jeffrey  B.  Burl,  Code  EC/B1  1 
Department  of  Electrical  &  Computer  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

6.  Prof.  George  Thaler,  Code  EC/Tr  1 
Department  of  Electrical  &  Computer  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

7.  Lt.  Kurtis  B.  Miller  2 
Navy  Undersea  Warfare  Center,  New  London  Detachment 

New  London,  Connecticut  06320 

8.  Harold  M.  Miller  1 
2314  Bennett  Ave. 

Glen  wood  Springs,  Colorado  81601 

9.  Kenneth  E.  Smith  1 
109  Bluegrass  Dr. 

Schererville,  Indiana  46375 


57 


DUDLEY  KNOX  LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY  CA    93943-5101 


GAYIORD  S 


